Constraints on Non-Gaussianity from Sunyaev— Zeldovich Cluster Surveys 



Daisy S. Y. Mak^'B and Elena Pierpaoli^^lll 

^Physics and Astronomy Department, University of Southern California, Los Angeles, California 90089-0484, USA 

We perform a Fisher matrix analysis to forecast the capability of ongoing and future Sunyaev- 
Zeldovich cluster surveys in constraining the deviations from Gaussian distribution of primordial 
density perturbations. We use the constraining power of the cluster number counts and clustering 
properties to forecast limits on the /nl parameter. The primordial non-Gaussianity effects on 
the mass function and halo bias are considered. We adopt self-calibration for the mass-observable 
scaling relation, and evaluate constraints for the SPT, Planck, CCAT-like, SPTPol and ACTPol 
surveys. We show that the scale-dependence of halo bias induced by the local NG provides strong 
constraints on /nl, while the results from number count are two orders of magnitude worse. When 
combining information from number counts and power spectrum, the Planck cluster catalog provides 
the tightest constraint with (t/j^l = ^ (68% C.L.) even for relatively conservative assumptions on 
the expected cluster yields and systematics. This value is a factor of 2 smaller than the la error 
as measured by WMAP CMB measurements, and comparable to what expected from Planck. We 
find that the results are mildly sensitive to the mass threshold of the surveys, but strongly depend 
on the survey coverage: a full-sky survey like Planck is more favorable because it can probe longer 
wavelengths modes which are most sensitive to NG effects. In addition, the constraints are largely 
insensitive to priors on nuisance parameters as they are mainly driven by the power spectrum probe 
which has a mild dependence on the mass-observable relations. 



I. INTRODUCTION 

Cosmological inflation has emerged as the most popu- 
lar scenario of the early universe and it predicts a near 
scale invariant power spectrum and close-to Gaussian 
distribution for the primordial curvature inhomogeneities 
that seeds large scale structure (LSS). Various inflation- 
ary models produce different levels of departures from 
Gaussianity. For example, the slow-roll and single field 
inflation model produces tiny amount of departures from 
Gaussianity, while other models predict sizable amount of 
primordial Non-Gaussianity (NG) that could be observed 
with current experiments such as CMB and galaxy clus- 
ters [e.g.[l|,[2]- Therefore, any detection of primordial NG 
would open a new and extremely informative window on 
the physics of inflation and the very early Universe. 

Current measurements of the CMB [e.g. Q and 
LSS [e.g. a found that the distribution of primordial 
fluctuations is consistent with Gaussianity, however, that 
bound is still several orders of magnitude away from test- 
ing primordial NG at the level predicted by slow-roll in- 
flation. Galaxy clusters are in principle sensitive and 
powerful tool for this purpose because they trace the 
rare, high mass tail of density perturbations. As a result, 
the changes in the shape and evolution of the mass func- 
tion of dark matter halos are most sensitive to departures 
from Gaussianity. The effect of NG on the mass function 
has been investigated in e.g. and was validated by 

large-scale cosmolo gica l simulations with non-Gaussian 
initial conditions [7l. llll - fl3| . More recently NG effects on 
the large scale clustering of collapsed haloes were studied 
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by [13, ■ They found that the linear biasing pa- 

rameter acquires a scale dependence, which modifies the 
power spectrum of the the distribution of cosmic struc- 
tures most prominently at large scales. This unique sig- 
nature serves as a powerful way to constrain and forecast 
the nature of NG assumption. For instance, predictions 
for the cluster abundance and clustering of galaxy clus- 
ters expected from various surveys with primordial NG 
conditions were presented in [13, • 

In this work, we aim to forecast the capability of on- 
going and future SZ cluster surveys in constraining pri- 
mordial NG, using the Fisher Matrix approach. Several 
initial studies have explored such possibility with future 
galaxy cluster surveys in different wavebands, e.g. X- 
ray [l^, [lO] and optical (2ll - [23l | , and they all showed that 
the constraints from galaxy clusters are quite strong. 
With the use of SZ cluster survey, we would complete 
the spectrum of cosmological applications with galaxy 
clusters. Forthcoming SZ experiments will provide large 
samples of mass selected clusters, and multi-frequency 
foUowup observations can add information on the clus- 
ter mass. In addition, these surveys will detect clusters 
at high redshift to test non-Gaussianity in the regimes 
where its effects are pronounced, that is, the high mass 
tail of the mass function and the large scale power spec- 
trum of the cluster distribution. This makes the SZ clus- 
ters more favorable since X-ray and optical clusters are 
not as efficient in detecting high redshift clusters. 

This paper is organized as follows. Wc begin by 
presenting the surveys and expected cluster samples in 
Sec.|lTl In Sec. IIIII we present the parametrization of pri- 
mordial non-Gaussianity effects on the halo abundance 
and clustering. Sec. IIVI details the Fisher formalism em- 
ployed here, as well as the fiducial cosmology adopted. 
The forecasted constraints are presented in Sec. |Vl We 
discuss our results in Sec. IVII and conclude in Sec. IVIIl 



A. The Planck Catalog 
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FIG. 1: Upper: Mass limit of cluster surveys. Lower: The 
redshift distribution of clusters in the Planck (black), ACTpol 
(blue), CCAT-like (red), SPT (green), and SPTpol (magenta) 
survey in the fiducial ACDM cosmology (solid). The bottom 
panel shows the fraction deviation from the gaussian number 
count when /nl = 100. 



II. CLUSTER SURVEYS 



Planck is imaging the whole sky with an unpreeedented 
eombination of sensitivity {AT/T ^ 2 x 10~^ per beam 
at 100 - 217 GHz), angular resolution (5' at 217 GHz), 
and frequency coverage (30 — 857 GHz). The SZ sig- 
nal is expected to be detected from a few thousand in- 
dividual galaxy clusters. Planck will produce a cluster 
sample with median redshift ^ 0.3 (see Fig. [1] upper 
left panel). The SZ observable is the integrated Comp- 
tonization parameter Y = J y dJlciustor out to a given 
radius. For Planck, a 5a detection threshold ensuring 
high level of completeness (about 90%) corresponds to 
^200, pc > 2 X 10~'^arcmin^ [2^|, where I200, pc is the inte- 
grated comptonization parameter within r2oo,p^, the ra- 
dius enclosing a mean density of 200 times the critical 
density. The early release from the Planck Collabora- 
tion gives a sample of 189 high signal-to-noise SZ clusters 
with > 6a detection. It is therefore likely that our as- 
sumed detection threshold will be eventually reached in 
future data releases. For an SZ survey, its flux limit can 
be translated into a limiting mass by using simulation- 
cahbrated scaling relations |29| : 
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In order to mitigate the effect of overestimation of 
unresolved clusters at low redshift, we further restrict 
^^iim,200pc to be at least 1O^^M0 at all z. With all 
these criteria, the Planck survey is expected to detect 
~ 1000 clusters. The mass threshold we find with this 
approach is consistent with the one in [33|- While we 
keep l200,pc = 2 x lO^'^arcmin^ as our reference mini- 
mum value for presentation of the main results, we will 
also discuss predictions for a lower mass threshold, corre- 
sponding to 1200, pc = 10~'^arcmin^. With such threshold, 
the completeness of the S/N > 5 sample is reduced to 
about 70% and the total number of clusters is 2700. 



We will investigate the predictions for the five surveys 
described below. While we try to obtain as realistic sur- 
vey specifications as possible, in particular for the mass 
limit as function of redshift Miiin(-z), the intrinsic scatter 
in the mass observation relations and the lack of previ- 
ous large samples of SZ clusters necessarily make these 
quantities somewhat uncertain. In particular, the rela- 
tion between cluster mass and SZ signal is still very un- 
certain (e.g. [23427I ). The final mass limits as a function 
of redshift are shown in the upper panel of Fig. (TJ and 
the resulting expected number of clusters for each survey 
is shown in the lower panel of Fig. [T] These limits are 
derived considering the characteristics of each survey, as 
specified in the following subsections. 



B. SPT and SPTpol 

The SPT survey is currently observing the sky with a 
sensitivity of 18^K/arcmin^ at 148 GHz, 218 GHz, and 
277 GHz. This survey covers w 2500 square degrees 
of the southern sky (between 20h > RA > 7h, -65° < 
6 < —30°) with a projected survey size and cluster mass 
limit well matched to the Stage HI survey specification 
of the Dark Energy Task Force 3]J. For the mass limits, 
we employ the calibrated selection function of the survey 
by [sij . This is based on simulations and used to provide 
a realistic measure of the SPT detection significance and 
mass. Disregarding the scatter in the fitting parameters 
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for this relation, we use here: 



D. CCAT-like 
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where ^ is the detection significance. For the SPT survey, 
we take clusters detected at ^ > 5 which ensure a 90% 
purity level. Currently, the SPT team is setting a low 
redshift cut at Zcut = 0.3 in their released cluster sample, 
due to difficulties in reliably distinguishing low-rcdshift 
clusters from CMB fluctuations in single frequency ob- 
servations. Nevertheless, with upcoming multi-frequency 
observations, a lower cut Zcut = 0.15 will likely be at- 
tained. We therefore apply this cut in our work. With 
this, the SPT survey is expected to detect ^ 500 clusters. 

In addition to this, we also consider the upcoming SPT 
polarization survey (hereafter SPTpol) which will have 
an increased sensitivity of 4.5/iK/arcmin^ at 150 GHz 
for a 3 year survey and sky coverage of 625 square de- 
grees. We scaled the mass limits by a factor of 3.01/5.95 
in Eq. ^ to match with the expected mass limits of 
SPTpol clusters (Benson 2011, private communication). 
We again use Zcut = 0.15, resulting in a total expected 
number of ~ 1000 clusters. While these are the limits we 
use for our main results, we also discuss outcomes that 
consider a lower mass limit, corresponding to ^ = 4.5 
(80% purity). With this mass limit, SPT would find 800 
clusters and SPTPol about 1400 clusters. 



C. ACTpol 

The Atacama Cosmology Telescope (ACT) has been 
observing a portion of the southern sky since 2008 con- 
sisting of two strips of the sky, each 4 degrees wide in dec- 
lination and 360 degrees around in right ascension, one 
strip is centered at (5 = — 5°, and the other is centered at 
5 = -55° nil. With a sensitivity of « 35^K/arcmin2, 
only about 100 clusters are expected to be detected. In- 
stead, we turn to the newly developing dual-frequency 
(150 GHz and 220 GHz) polarization sensitive receiver 
(hereafter ACTpol [l^] and reference therein) to be de- 
ployed on ACT in 2013. One of the three ACTpol observ- 
ing seasons will have a wide survey covering 4000deg^ to a 
target sensitivity of 20/.iK/arcmin^ in temperature at 150 
GHz. With the wide field, they aim to find ~ 600 clusters 
in the ACTpol survey. The survey is 90% complete above 
a limiting mass of Miim.200p = 5 x lQ^'^MQh~^ (Sehgal 
2011, private communication), and we therefore assume 
this as our redshift-independent mass limit for ACTpol. 
As in SPT, the ACT team also put a low redshift cut 
in their parameter determination works and we likewise 
take Zcut = 0.15 for ACTpol, resulting in a total expected 
number of ^ 500 clusters. We also present in the discus- 
sion section the results corresponding to a lower mass 
limit; Miim.200p = 4 x 10^'* Mq/i"!, which would result in 
a catalog of about 1000 clusters. 



The CorneU Caltech Atacama Telescope (CCAT-like) 
will be a 25 meter telescope for observations at submm 
wavelengths. It will combine high spatial resolution 
(w 5"), a wide field of view (20'), and a broad wavelength 
range (100 GHz-405 GHz) to provide an unprecedented 
capability for deep, large area, multi-color SZ surveys of 
galaxy clusters to complement narrow field, high resolu- 
tion studies with ALMA. The telescope aims for initial 
observations in 2015 and will be the first to provide a 
large sample of clusters with high spatial resolution SZ 
profiles that aid in studying cluster astrophysics. 

We consider here a prototype of CCAT-like catalog 
from the CCAT-like LWcam survey which will cover 2000 
square degrees of the sky for 1000 hours, resulting in a 
noise level of 12^/ir/arcmin at 220 GHz. Following the 
LWcam Design Study Proposal (private communication, 
Vandcrlindc 2012), wc estimate its mass limits using the 
framework developed for the SPT (see Sec. Ill B]) . Like- 
wise, we consider clusters detected at ^ > 5 (90% purity 
level) in the redshift range z = — 1. Using Eq. (I2|), the 
CCAT-hke catalog is expected to detect 1500 clusters. 



III. PRIMORDIAL NON-GAUSSANITY 

Early universe models predict deviation from Gaussian 
initial conditions. Primordial non-Gaussianity induced 
by inflationary models can be conveniently parametrized 
by a nonlinear coupling parameter /nl and the Bardeen's 
gauge invariant potential $. This can be written as the 
sum of a linear Gaussian term and a non-linear second 
order term that encapsulates the deviation from Gaus- 
sianity [33l - [36j 



* = $g + ./nl(*g-(*g)) 



(3) 



The parameter /nl determines the amplitude of the 
non-Gaussianity. In this work, we adopt the large scale 
structure convention for defining the fundamental param- 
eter /nl • As noted by [ll| , the primordial value of $ has 
to be linearly extrapolated at z = in the LSS con- 
vention, so that /nl = 9{z)lg{0)f§f^'' ~ ^■SfSIf'' for 
the ACDM model. We note that the factor 1.3 is ap- 
proximate since the growth factor g{z) is cosmological 
dependent. This means that any constraints gathered 
from CMB data should be increased by 30% in order to 
comply with the convention adopted here. 

We define here several terms which appear in subse- 
quent text. The relation between the power spectrum 
matter density fiuctuation extrapolated at z = 0, P{k), 
and the power spectrum of the Newtonian potential, 
P$(fc) is 
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P$(fc) = A4(fc)P<t,(fc) (4) 
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The primordial matter power spectrum is scale free, i.e. 
P{k) = Ak"", and therefore the potential power spectrum 
can be written as P^{k) = £i^|^A:»-4 = Bfc»-4, 

In the case of non-Gaussianity, the random field of the 
potential $ cannot be described by the power spectrum 
P$ — Bk"^^ alone. Higher order moments, in particular 
the bispcctrum _B$(fci, k2,kz), are required. The bispec- 
trum is defined on the basis of the Fourier transform of 
the three point correlation function (|$(fci)$(fc2)$(A:3 
as follows 



$(fcl)$(fc2)$(fc3); = (27r)^(5D(fcl + fc2 + A:3)S$(fci,fc2,fc3) 

(5) 

The shape of the non-Gaussian bispectrum is related 
to the fundamental physics of the early universe and the 
evolution of the inflation field. A wide class of inflation- 
ary scenarios lead to non-Gaussianity of the local type 
in which the bispectrum of the Bardeen's potential is 
maximized for squeezed configurations. Example of this 
scenario is the curvaton model which involves an addi- 
tional contribution to the curvature perturbations by a 
light field [13 . The parameter /nl of the local type is a 
constant in space and time with /nl <C 1, and is expected 
to be of the same order of the slow-roll parameters [s^ . 
In such case, the bispectrum has a simple form [s^: 
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for the Gaussian mass function na^M, z). Using LMVJ 
approach, the correction factor R is 



i?(A/,z) = l+l^ 
6 dc 



54 ^2 
^3(^-2^-1) 



^^(4-1) 

d In (TM (Jm 
(8) 



where 5c = Ac/D{z) is the critical density for collapse, 
D{z) is the linear growth factor, ctm is the rms of pri- 
mordial density fluctuations on the scale corresponding 
to mass M, S3 = /nlM3(-^-'^)/o'm t^^*^ normalized skew- 
ness, ^3 is the third order moment. For the local non- 
Gaussianity, /Z3 can be computed as 
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where kf2 = k\ + k"^ + 2fikik2- This integral is compu- 
tational intensive. In order to reduce the workload of 
this calculation, we instead use the fitting formula of S3 
by jisj which is shown to have sub-percent accuracy: 
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We verified the accuracy of this fitting formula by di- 
rectly comparin g th is with the numerical table publicly 
available online [53, and we found that they agree with 
each other to percent level. 



Mass Function 



B. Halo Bias 



Non-gaussianity lead to a modified mass function with 
respect to the gaussian case and are usually expressed 
as perturbations of the gaussian mass function. There 
are several prescriptions of the corrections in mass func- 
tions of collapsed objects (e.g. In this work, we 
adopt the approach of [40] (hereafter LMVJ), in which 
the probability distribution for the smoothed dark mat- 
ter density field is approximated using the Edgeworth ex- 
pansion truncated at the the first few orders. The LMVJ 
approach was shown to give reasonable agreement with 
full numerical simulations of structure formation (Tl| . 
provided that the linear overdensity threshold for col- 
lapse is corrected for ellipsoidal density perturbations, 
i.e. Ac — > ^cy/q where q « 0.75. 

In this prescription, the non-Gaussian mass function 
n-MciM, z) can be written as a function of a Gaussian 
one, TT-G, multiplied by a non Gaussian correction factor 
R{M,z), 



mciM, z) = R{M, z)nG{M, z) 



(7) 



where R{M,z) = nG,ps/"NG,PS, and riQ^ps and uncfs 
are the Gaussian and non-Gaussian mass function respec- 
tively computed according to the Press and Schechter for- 
mula [4l| . In this work, we adopt the formula from 42 1 



The halo bias acquires an extra scale dependence due 
to primordial non-Gaussianity of the local type [l^, i.e. 



6ng(M, z, k) = 6g(M, z) + Ab{M, z, k) (11) 



where 



Afo(Af,z,fc) = [6g(M,z) - l],5,r;^(fc) 



(12) 



The term r^(fc) encapsulates the dependence on the 
scale and mass. For the local non-Gaussianity it can be 
written as 



r^^Xfc) = ^^^^i^J^ dk,kl^Uk^)p{k^) 



^\ d^MR{ki2) 



P(fcl2) 

P{k) 



(13) 



where fc^j = kf + k^ + 2^kik. The function Tpi{k) is fiat 
at small scales but scales as k~^ at large scales {k > 0.01 
Mpc/h), so that a substantial deviation in the halo bias is 
expected at those scales. We adopt the formula from [i^l 
for the Gaussian bias Bq. 
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IV. FISHER MATRIX FORECAST 



The Fisher information Matrix (FM hereafter) is de- 
fined as 



q/3 
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where £ is the hkehhood of a data set, e.g. a cluster sam- 
ple, written as a function of the parameters Pa describing 
the model. The parameters pa comprise the cosmologi- 
cal model parameters as well as "nuisance" parameters 
related to the data set (e.g., mass calibration). 



A. Cosmological parameters 

Throughout this paper, we assume a spatially flat 
(ilfc = 0) cosmology. Our model comprises a total of 
seven cosmological parameters and the non-gaussianity 
/nl parameter which are left free to vary. The seven pa- 
rameters and their fiducial values (in parenthesis, taken 
from the best-fit flat ACDM model from WMAP 7yr 
data, BAO and Hq measurements Q) are: baryon den- 
sity parameter fib/i^ (0.0245); matter density parameter 
uj„i = ^m.h^ (0.143); dark energy density Q,\ = 1 — ilm 
(0.73); power spectrum normalization erg (0.809); index 
of power spectrum Us (0.963); effective dark energy equa- 
tion of state through w{z) ~ wq + {1~ a)wa, with fiducial 
values wq = —1 and Wa — 0- The Hubble parameter is 
then a derived parameter given hy h = •\/wm/(l — ^^a) — 
0.73 in the fiducial case. 

In the following, we first discuss the Fisher matrix 
for number counts and clustering of clusters, before de- 
scribing the calibration parameters and CMB priors. 
Throughout, we divide the redshift range into bins I of 
width Az = 0.05. Further, we bin clusters in logarithmic 
mass bins m of width AlnM = 0.3 from the minimum 
mass M\i„i{z) for each survey (Sec.lll]) up to a large cut- 
off mass of Minax = lO^^A^Q. Since the mass limit varies 
with redshift, the number of mass bins thus also varies 
somewhat across the redshift range. 
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where AQ is the solid angle covered by the cluster sur- 
vey, In Mi,rn = In ^-^11111(2^;) + mA In M, and n{M, z) is the 
mass function given in Eq. ([7]). Following Lima and Hu 
[45| . we take into account the intrinsic scatter in the re- 
lation between true and observed mass, as inferred from 
a given mass proxy, by the factor p{M°^\M) which is 
the probability for a given cluster mass with M of hav- 
ing an observed mass . Under the assumption of 
a log-normal distribution for the intrinsic scatter, with 
variance cf^ ^ , the probability is 
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With these notations, we parameterize the M°^ — M rela- 
tion, in addition to the intrinsic scatter, by a systematic 
fractional mass bias Bm- With this prescription, the final 
expression for the number count FM is: 



Af2Az d'^V 



(19) 
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where erfc(x) is the complementary error function. 



C. Power spectrum 



B. Number counts 

The FM for the number of clusters Ni,„i within the 
Z-th redshift bin and m-th mass bin is 



dNi.ra dN, 
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where the sum over I and m runs over intervals in the 
whole redshift range z = — 1 and cluster mass range 
[Mii^{z),Qo]. We can write the abundance of clusters 
expected in a survey, within a given redshift and mass 
interval, using the mass function as: 



We define the FM for the power spectrum of galaxy 
clusters as 
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where the sum over m,n runs over mass bins, while 
the sum in I and i runs over z and k intervals respec- 
tively. In what follows, we use fcmax = 0.1 /iMpc~^ 
with Alogj^gfc = 0.017 for all surveys. Primordial NG 
of the local type modifies the shape of the power spec- 
trum of galaxy clusters by introducing a scale-dependent 
bias on very large scales, therefore the choice of largest 
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scale that can be probed by a cluster survey would signif- 
icantly affect the constraints on /nl- In practice, this is 
taken into account by introducing the window function 
of the survey in the calculation of P{k). However, for 
simplicity, we approximate the influence of the window 
function by considering a cut off in the choice of fcmin- 
We use /cmin = 10~^ for the all-sky Planck survey, and 
fcmin = 10~^ for the partial sky surveys. 

P^"{ki,zi) is the cluster cross-power spectrum for 
mass bins m and n, calculated for the given redshift and 
wavenumber through 



(21) 



Here, is the mass function weighted effective bias. 
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(22) 

The effective volume for mass bins m, n, wave number 
ki, and redshift zi is given by (see App. A of [46j ) 
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where Vo{z) is the comoving volume of the redshift slice 
[zi — 0.01, -I- 0.01] covered by the given survey, and 
nmizi) is the cluster number density for mass bin m at 
redshift zi. The effective volume gives the weight car- 
ried by each bin in the (z, k) space to the power spec- 
trum Fisher matrix, and hence quantifies the amount 
of information contained in a given redshift- and fc-bin. 
Fig- m shows the redshift and scale dependence of the ef- 
fective volume for the five cluster surveys. We find that 
Vcff < 0.3 Vq for all redshifts and surveys considered, even 
when not binning in mass, hence the cluster power spec- 
trum is shot-noise dominated for all surveys. As the lower 
panel of Fig. [2] illustrates, Planck is most limited by shot 
noise, while SPTpol is least limited, as expected from 
their respective mass limits and coverage. 



D. Calibration parameters 

In self-calibrating the true and observed cluster mass 
(Eq. (|18|) ). we introduce four nuisance parameters which 
specify the magnitude and redshift-dependence of the 
fractional mass bias Bm(z) and the intrinsic scatter 
<Jin m{z). Following [4^|, we assume the following 
parametrization: 



BM{z) = BMo{l + zr 
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FIG. 2: The dependence on redshift (top) and wavenumber 
(bottom) of the effective volume (Eq. (|23l) ) for a single mass 
bin and each survey: Planck (black), SPT (green), SPTpol 
(magenta), ACTpol (blue), and CCAT-like (red). The effec- 
tive volume is a weak function of wavenumber k but strongly 
depends on the redshift. 



Therefore the four nuisance parameters are Bmq, ol, 
<y\nM,o, and /3. A negative value for Bm corresponds to 
an underestimation of mass. The mass bias accounts for 
the possibility of a systematic offset in the calibration of 
the observable mass scaling relation. We adopt fiducial 
values of i?j\/o = 0, a = 0, (7\nM = 0.1, = 0. In deriv- 
ing the main results, we will not make any assumption 
on the four nuisance parameters and leave them free to 
vary. We will study the effect of assuming different priors 
on the four nuisance parameters on the /nl constraints 
in Sec. lyiAl 



E. CMB Prior 

In the following, we present results with the Fisher 
matrix for the Planck CMB temperature power spec- 
trum Ci added to the constraints from cluster counts 
and power spectrum. We calculate the full CMB fisher 
matrix with CAMB [13 and method described in [isi. 




For the Planck experiment, we use the three frequency 
bands 100, 143 and 217 GHz, and the C; are calcu- 
lated up to ^max = 2500. Our fiducial parameter set for 
the CMB experiment is, as described in the DETF re- 
port nil, e = {ns,^bh'^,^A,^nih'^,wo,As,T), where As 
is the primordial amplitude of scalar perturbations and r 
is the optical depth due to reionization. After marginal- 
izing over the optical depth, we transform the Planck 
CMB fisher matrix to our cluster survey parameter set 
^' = {n's,ftbh^,^A,^mh'^jWo,as) by using the appropri- 
ate Jacobian matrix. The CMB imposes strong prior 
on the cosmological parameters. For example, ftm.h^ is 
known to be measured with the CMB power spectrum to 
an exquisite precision, and this helps in breaking param- 
eter degeneracies in the constraints from cluster surveys. 
On the other hand, we note that the CMB power spec- 
trum docs not add constraints on /nl and therefore we 
compute the CMB fisher matrix for the Gaussian pertur- 
bation. 



V. RESULTS 




0.62 r 




0.0100 0.0148 0.0219 0.0324 0.0479 0.0703 




0.20 0.64 1.08 1.96 7.40 ?.83 



FIG. 3: Relative deviations of tlie non-gaussian power spec- 
trum from the gaussian power spectrum, i.e. AP/Pq for the 
Planck survey and /nl = 100. 



A. Cluster counts and Power Spectrum 

The first four rows of Tab. U summarize the marginal- 
ized /nl constraints from dN/dz. Under the standard 
setup: assuming we have spectroscopic redshifts Az = 
0.05 and binning in cluster mass of AlogAf = 0.3, the 
/nl constraints are ~ 10"^. The relative constraining 
power of the different surveys can easily be interpreted 
by looking at AA^ shown in Figure [TJ The Planck and 
CCAT-like survey, which have the largest AA^ at z < 0.3 
and z > 0.3 respectively, similarly give the tightest con- 
straints of (T/pjL = a,nd CT/pjL = 815 respectively at the 
68% CL. The other three surveys are less stringent and 
the constraints are > 80% worse than CCAT-like. 

We study how the constraints change if we ignore the 
redshifts and consider combining information on all red- 
shifts. With only one redshift bin, none of the surveys 
is able to constrain /nl- This implies redshift informa- 
tion on each clusters is essential in constraining /nl using 
number counts. So far we have assumed the optimistic 
scenario in which the cluster redshifts are spectroscopic. 
We relax this choice by using the photometric redshifts 
(i.e. Az = 0.1), this worsens the constraints, particularly 
for Planck in which cr/p^L increased by 59% . 

The constraints from power spectrum only arc summa- 
rized in the row 5-9 in Tab. [D Overall the constraints are 
almost two orders of magnitude better than that from 
number counts only. This is because the scale depen- 
dence of the halo bias is very sensitive to non-gaussianity 
(e.g. ti2.) JLi, ,.5^ ) . The best constraint come from the all- 
sky Planck survey, with (T/„l = 7. This result is largely 
due to its large sky cov erage ( .fskv ~ 0.7). It was shown 
in previous studies [e.g. [Tl|, |40] that the deviation of the 
halo bias, and hence the power spectrum, between non- 
gaussian and gaussian case is most prominent at large 



scale. We illustrate this in Fig. [3] with the AP/P in 
different redshift and wavelength bins. The survey that 
probes longer wavelength modes, as in the case of Planck 
with kmin = 10~^, therefore has the highest sensitivity to 
/nl- 

We here quantify the sensitivity of the non-gaussianity 
constraints to the adopted fcmin value in the analysis of 
the power spectrum. We address this by computing a/^,^ 
as a function of /cmin, as shown in Fig. HI As expected, 
the constraints significantly improve when we consider 
smaller k values. This is partly because more informa- 
tion is included in the fisher matrix and, more impor- 
tantly, the effect of NG on halo bias is most prominent 
at the largest scales. It is interesting to note that if we 
limit fcmin to 10~^ Mpc/h also for the Planck survey, then 
the derived /nl constraint would be similar to, and even 
slightly worse than the one derived from the CCAT-like 
and SPTpol survey. 

It is worthwhile to note that the effective volume Kff 
also impacts the /nl constraints by its z and fc depen- 
dence on the fisher matrix, as shown in Fig. [5] This is 
particularly useful in understanding the relative merits of 
the four partial sky surveys that have similar survey area 
and fc-range. Indeed the relative constraining power of 
these four surveys is reflected in the redshift distributions 
of the effective volume: while CCAT-like and SPTpol has 
the largest Vcs at z < 0.7 and z > 0.7 respectively (upper 
panel of Fig. [2]), their errors on /nl are similar and are 
40% smaller than ACTpol and SPT. To further under- 
stand which redshift range contributes most to the /nl 
constraints, we show in Fig. |4] the /nl constraints as a 
function of the maximum cluster redshifts. If we limit the 
cluster samples to z < 0.6, both CCAT-like and ACTpol 
survey give tighter constraint than SPTpol and SPT. As 
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TABLE I: Marginalized la errors on /nl- The labels in the second column means the following, standard: The standard 
setup as indicated in Sec. lIVBl and Sec. IIV Cl less conservative: The conservative mass limit of the cluster survey is considered 
(see Sec. IV Bl photo-z: Redshift binning correspond to the error of photometric redshift, i.e. Az = 0.1; one mass bin: No mass 
slicing in the fisher matrix, i.e. only one mass bin. 



Probes Planck ACTpol SPT SPTpol CCAT-like 



dN/dz 


standard 


987 


1473 


2319 


1463 


815 




less conservative 


497 


1069 


1810 


1230 


- 




photo-z 


1540 


1588 


2355 


1462 


912 




one mass bin 


CXD 


CO 


CXD 


oo 


oo 


P(k) 


standard 


7 


24 


21 


15 


13 




less conservative 


4 


16 


17 


13 






photo-z 


8 


24 


20 


15 


13 




one mass bin 


160 


272 


267 


124 


115 


dN/dz 


+ P(k) standard 


7 


24 


20 


15 


13 




less conservative 


4 


16 


16 


13 






photo-z 


7 


24 


20 


15 


13 




one mass bin 


156 


269 


255 


119 


113 



we extend to larger maximum redshift, SPTpol survey 
then gives tighter constraint than SPT and ACTpol sur- 
vey as its Vee{z) is decreasing at a smaller rate. This 
result shows that the constraints leverage on clusters in 
disjoint redshift ranges and these surveys provide com- 
plementary information on /nl constraints from power 
spectrum. 



B. Combined constraints 

The constraints on /nl when combining both number 
counts and power spectrum are summarized in rows 10-13 
of Tab. m While combining the two probes helps to break 
degeneracies and improves constraints on cosmological 
parameters, the constraint on /nl does not improve. The 
constraining power is mainly driven by the power spec- 
trum. We note that the best constraint estimated here 
(^/nl ^ from Planck cluster survey) is a factor of 2 
(when accounting the difference in f^f^ and f§ff^) bet- 
ter than the current upper limits as measured by CMB 
experiments (e.g. 32 ± 21 and a factor of 14 better 
than large scale structure probes (e.g. A/nl = ±96 Q). 
In the long run, however, constraints from new CMB 
measurements from Planck would push the upper lim- 
its further down to a/^^ < 5 [5l|. These Planck limits 
also consider polarization maps and disregarding poten- 
tial foregrounds, thus it these CMB limits should be re- 
garded as the optimistic constraint. Next generation all- 
sky X-ray surveys ^19.. 2{)\wouldachieve<7fj.-,^ w 10. Note 
that constraints from CMB are based on different physics 
than large scale structures and suffered from different 
systematics. The SZ clusters surveys therefore provide a 
useful complement to the /nl information we can derive 
from the CMB maps. 




0.001 0.010 0.100 

k^,„ (Mpc/h) 




0.0 0.2 0.4 0.6 0.8 1.0 

Zmax 



FIG. 4: Fully marginalized constraints on /nl from the power 
spectrum of clusters only, as a function of minimum wavenum- 
ber fcmin (upper) and maximum cluster redshift Zmax. (lower). 



Fig. [U illustrates the most important degeneracies of 
/nl with cosmological parameters (erg, Wq, Wa) for the 
Planck survey. These plots clearly demonstrate the com- 



9 



plementarity that the number counts and clustering have 
to constrain /nl and standard cosmological parameters, 
particularly the redshift dependence part of equation of 
state of dark energy Wa- It is obvious that the /nl is 
almost non-degenerate with other cosmological parame- 
ters. Only mild degeneracies are seen in the joint con- 
straint with wq and Wa- Therefore the constraints on 
other parameters have little effect on /nl- 

It is important to keep in mind that these results are 
based on a conservative mass limits, i.e. clusters are ex- 
pected to be detected with S/N > 5 for all surveys. We 
examine improvements in the /nl constraints when using 
more optimistic mass limits for each survey according to 
what is outlined in SecHIl Using number counts only, the 
constraint from Planck is improved by 50%, while con- 
straints from other surveys are only slightly improved 
(by 16 — 27%). However, these constraints are worse 
than those from power spectrum. Using power spectrum 
only, the constraints are improved only marginally for 
SPTpol and SPT, but largely by 30 - 40% for Planck 
and ACTpol. It is interesting to notice that ACTpol is 
now slightly better in constraining /nl than SPT while it 
was the opposite case when using the conservative mass 
limits. This is mainly due to the larger increase in de- 
tected clusters for the ACTpol survey in the optimistic 
mass limits. 

As a worse case scenario, we also consider the case 
when there is no information on the cluster mass, i.e. 
one mass bin. For all surveys, the combined constraints 
are significantly worsened by a factor of > 10. Yet, the 
CCAT-like and SPTpol constraints with single mass bins 
arc still comparable to current upper limits. 



VI. DISCUSSION 

A. Uncertainties in scatter of mass observable 
relations 




).a04 0.806 0.808 0.81 



0.812 0.814 



dN/dz + Planck 

E(k) + Planck 

+ Planck ; 




-0.6 -0.4 



dN/dz + Plane < 
_Em,t£lans< 




-6 -4 -2 2 4 6 



Our forecasts of /nl assumed no prior on the nuisance 
parameters in the mass-observable relations: Bm,o, cj\/,Oj 
a, /?. However, allowing more freedom to the scaling 
relations parameters results in degradation of the final 
constraint. In practice, we can expect some external 
constraints on these parameters by using detailed stud- 
ies of individual clusters or combining different informa- 
tion from optical, lensing. X-ray and SZ measurements. 
To better understand the importance of self-calibration 
of systematic uncertainties on the non-Gaussianity con- 
straints, we repeat the forecasts on different priors on the 
four nuisance parameters. 

The results are summarized in Tab. HH We have 
used the current knowledge on the calibration on the 
mass proxies from X-ray and lensing measurements [e.g. 
[5^ and assume priors on the nuisance parameters of 
Actm.o = 0.1, A/3 = 1, = 0.05, and Aa = 1. 

We refer this to the weak prior. Furthermore, we refer 
to the case strong prior when we assume that all these 



FIG. 5: Joint constraints on the /nl and (counterclockwise 
from top left) as, Wo, and Wa- All curves denote 68% con- 
fidence level, and are for number counts only (blue), power 
spectrum only (cyan), and combination of the two (green). 
The Planck CMB power spectrum priors are assumed. 



four parameters are held fixed at their fiducial values. 
We find that the /nl constraints only slightly improves 
when using number counts only and has negligible im- 
provement in the joint constraints when considering the 
weak prior. In the case of strong prior, the constraints 
are significantly tightened by a factor of 4 — 8 from num- 
ber count only, but negligibly from power spectrum only. 
This results in appreciable improvements in the joint con- 
straints for the partial sky surveys but not for Planck, as 
measured from the cluster. 

As the /nl joint constraint is essentially driven by the 
power spectrum, we conclude that it is largely insensitive 
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to priors on nuisance parameters. This is consistent with 
the scenario shown in [l^ as they calculated the ratio of 
non-Gaussian to Gaussian at a wide range of values of 
the mass bias i?j\/o and scatter a mo- They showed that 
the ratio on number counts change by less than 0.1%, 
and is negligible on the effective bias except at the very 
large scales (40% at k fa 10~^Mpc~^). This justifies the 
weak dependence of the /nl constraints on prior since the 
choice of prior affects more the prediction for the cluster 
counts on /nlj thus show smaller dependence than that 
from clustering measurements. 



mass bins. This makes our constraint comparable to 
theirs even if our samples have much lower statistics 
((T/j^L = 7 for w 1000 clusters in the Planck survey vs 



11 for - 10^ clusters in the WFTX Wide sur- 



vey). On the other hand, if we do not consider mass 
bins in P{k) then our constraint degrade to (7/^^ = 156, 
which is an order of magnitude worse. The improvement 
mainly comes from the mass dependence of the halo bias 
which in turn provides extra information of the shape of 
the power spectrum. This suggests the slicing in mass 
bins compensates for the poor statistics. 



B. Comparison to previous work 



VII. CONCLUSIONS 



Forecast of the local-type non-gaussianity was previ- 
ously done with galaxy cluster probes in other wave- 
lengths, e.g. future X-ray surveys ([l^ and [lOl) and 
optical surveys (U, 22 1, [l^l). This study is the first 
attempt in exploring the /nl constraints from the SZ 
cluster survey and our results are in broad agreement 
with these previous studies. A detailed comparison be- 
tween our work and these studies is not straightforward 
because of the different survey specifications, prescrip- 
tions for the mass function and halo bias, and the type 
of cluster probes used. The more interesting compari- 
son can be made with the results from [l^ and [l^l who 
put constraints with cluster number counts and cluster- 
ing and similar subsamples. 

In terms of cluster sample ever considered by other 
authors, the subsample from cRosita [l^] of the most 
massive 1000 clusters (Mgoccrit > 2.2 x IO^^^/i-^Mq at 
z > 1, "magnificent 1000") is the most similar to the 
Planck cluster sample. This is because the limiting mass 
of the Planck catalog behaves in a more similar way to 
the X-ray catalog due to its large beams that signifi- 
cantly smoothes the signal, especially when the angular 
size of the cluster is small. We forecast a tighter con- 
straint (cT/^L = 7) than theirs (ct/^l = 26), when both 
assuming photometric redshift information and adding 
the Planck CMB prior. One should keep in mind that 
the non-gaussianity effect is more prominent at high red- 
shift, so the "magnificent 1000" which contains cluster 
at z > 1 should in principle be more sensitive to /nl 
constraint. The discrepancies could probably be due to 
a number of differences in the analysis: (a) they consid- 
ered different types and number of nuisance parameters 
that used to model the uncertainties in the scaling rela- 
tions, such as Lx — M, Tx — M. This is a conservative 
approach but at the same time degrade their constraints 
considerably, (b) the mass-redshift distribution of the 
X-ray sample are different from our SZ samples. The SZ 
samples tend to have higher fraction of massive clusters, 
i.e. the regime where the NG effects are more promi- 
nent. This makes the SZ clusters more sensitive to the 
/nl constraints. 

In terms of cluster probes, we use the same ones in jToj . 
except that we additionally consider slicing the P{k) in 



In this work, we exploited the large cluster samples ex- 
pected from current and upcoming SZ surveys to place 
constraints on primordial non-Gaussianity of the local 
type. Making use of the cluster number counts and power 
spectrum, and taking into account the self-calibration 
of mass-observable scaling relations, we employed the 
Fisher matrix analysis to forecast the sensitivities of var- 
ious SZ surveys in constraining the /nl parameter. 

The main results are presented in Sec. |Vl We find 
that the induced scale-dependence of halo bias (through 
a term that is proportional to fc~^) by local type non- 
Gaussianity provides a very effective way to put strong 
constraints on /nl- This makes the power spectrum a 
more powerful probe than number counts, by at least 
two orders of magnitudes. As a result, all-sky surveys 
such as Planck, which probe better the most NG sensi- 
tive regime at large scales, are more favorable to measure 
/nl- The best constraint we obtain is from Planck (com- 
bined constraints, (Jf^^ — 7) and is mainly driven by the 
power spectrum. The partial-sky surveys SPT, SPTpol, 
ACTpol, and CCAT-like, however, are a factor of > 2 
less constraining than Planck with (7/^^^ = 13 — 24. The 
best constraint we obtained, which is based on conserva- 
tive assumptions on uncertainties in mass-observable re- 
lations and mass thresholds, is a factor 2 better than that 
measured from WMAP CMB and is comparable to the 
expected results from the Planck CMB non-Gaussianity 
studies (/nl < 5). The SZ cluster surveys therefore pro- 
vide a useful complement to the /nl information we can 
derive from CMB maps. However one should note that 
these results are based on different physics and suffering 
from different systematics than the probes considered in 
this work. Wc also show that /nl have little degeneracy 
with other cosmological parameters and it is only mildly 
degenerate with wq and Wa- Thus the constraints on 
other parameters have little effect on /nl- 

We investigate the sensitivity of our results to various 
aspects of survey specification in Sec.lVIl We find that the 
errors on /nl are mainly driven by the power spectrum 
and therefore they are insensitive to priors on nuisance 
parameters. Only when we have perfect knowledge on 
the uncertainties of these parameter we can improve the 
constraint by a factor of 1.7-3.5, with Planck being the 
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TABLE II: Fractional improvement , with various priors (see Sec. 



Probes prior Planck ACTpol SPT SPTpol CCAT-like 



dN/dz 


weak 


1.63 


1.54 


2.13 


2.22 


1.63 




strong 


7.68 


3.48 


5.41 


4.39 


4.39 


P(k) 


weak 


1.00 


1.01 


1.01 


1.01 


1.01 




strong 


1.02 


1.04 


1.01 


1.02 


1.04 


dN/dz+P(k) weak 


1.01 


1.02 


1.01 


1.01 


1.01 




strong 


1.66 


3.46 


3.15 


2.80 


3.03 



least benefited. This, however, requires very good under- 
standing of tlie mass calibration which may be possible 
with multi-frequency foUowup measurements. In addi- 
tion, the cluster selection function plays a less significant 
role in which 0"/;.^^ is marginally improved when we con- 
sider more optimistic mass limits. 

Our results are in broad agreement with previous stud- 
ies that explores the constraining power of /nl from X~ 
ray and optical cluster surveys. The values of erf j^i^ we ob- 
tained are comparable to the ones quoted in these other 
works even if the SZ samples contain at least two orders 
of magnitude less clusters. This is because, unlike pre- 
vious authors, we considered mass slicing in power spec- 
trum when calculating the fisher matrix and this greatly 
improves the constraints. This suggests the slicing in 
mass bins compensates the poor statistics of the SZ clus- 



ter samples. We realize that the SZ surveys we consid- 
ered here are very promising and their cluster catalogues 
would be available in the very near future, since most 
of them are operating. Therefore, our /nl forecasts can 
readily be referenced. 
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